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Abstract 

Ongoing glacial retreat is expected to lead to numerous changes in glacier-fed 
rivers. This study documents the development of community composition of 
the hitherto widely neglected micro- and meiobenthos (MMB: bacteria, fungi, 
algae, protists, and meiofauna) in glacier rivers in response to the distinct habi- 
tat conditions driven by different stages of (de)glacierization. Our model is 
based on the glacier catchments of the Moll River (MC) and Kleinelendbach 
stream (KC), in the Austrian Alps, with 60% and 25% glacierization and glacier 
retreats of 403 and 26 m, respectively, since 1998. Analyses of overall catchment 
diversity and resemblance patterns showed that neither intense glacierization 
nor rapid deglacierization were predominant MMB determinants. This was 
ascribed to the specific environmental conditions at the MC, where the rapidly 
retreating Pasterze glacier has formed a harsh unstable proglacial, but also a 
benign floodplain area, with the former suppressing and the latter supporting 
the structural development of the MMB. Comparisons of similarly aged riverine 
habitats of the MC proglacial and the KC main channel further evidenced 
developmental suppression of the MMB (64 taxa) by the rapidly retreating MC 
glacier, unlike the moderate glacial retreat in the KC (130 taxa). Habitat condi- 
tions interacting with melt periods explained the differences in MMB resem- 
blance patterns, which themselves differentially reflected the spatiotemporal 
habitat settings imposed by the different glacier activities. The varying glacial 
influences were represented by a glaciality index (GIm) based on water temper- 
ature, electrical conductivity, and stream bed stability. The taxonomic richness 
of nematodes, rotifers, algae, and diatoms was distinctly related to this index, as 
were most MMB abundances. However, the strongest relationships to the GIm 
were those of nematode abundances and maturity. Our observations highlight 
the intense response of the MMB to ongoing glacier retreat and the utility of a 
simple index to reveal such patterns. 



Introduction 

Riverine habitats of glacier catchments are among the 
many highly vulnerable habitats that are threatened by 
present-day glacier retreat patterns and global climate 
change (Woodward et al. 2010; Vanham 2012). This 
threat becomes even more alarming given the far reaching 
consequences expected for riverine biota, whose complex- 
ity and complex relationships still constitute major con- 
temporary "black boxes" in ecological global change 



research. In general, river research has suffered from an 
imbalance in that it has largely focused on just a few 
organismal groups, mostly belonging to the macrofauna 
(Milner et al. 2001, 2009; Brown et al. 2009, 2010; Finn 
et al. 2010; Parker and Huryn 2011; Jacobsen and Dangles 
2012), rather than on complex community patterns (Sten- 
dera et al. 2012). However, baseline data for all glacier 
river inhabitants are required, not only to ensure reliable 
evaluations of the potential response of riverine biota to 
glacial and climate changes but also to enable progressive 
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and therein better targeted research into the resulting 
altered relationships between biota components. This in 
turn will allow estimations of ecosystem performances, 
now and in the future (Lecerf and Richardson 2010). 

Glaciers act as a major determinant of their catch- 
ment's riverine habitats, where a large degree of glacier- 
ization (Fiireder 2007) as well as ablation activity 
accounts for the harshness within these environments 
(McGregor et al. 1995) While glaciers and their climate 
responses have become better understood over the course 
of decades of research, which predicts ongoing glacier 
retreat (WGMS 2008; Dessler and Parsons 2010), uncer- 
tainties remain due to the complex reciprocity between 
climate, glaciers, and their surroundings (Pelto and 
Brown 2012). The uniqueness of each glacier is further 
exemplified by the variety of developmental stages, hydro- 
logical regimes, and habitat heterogeneity (templates and 
patchiness) within alpine riverine networks following 
glacial retreat in response to ongoing climate change. In 
fact, major changes are expected for catchments with a 
high degree of glacierization and for glaciers with signifi- 
cant ablation zones (Hook 2005; Dyurgerov et al. 2009; 
Farinotti et al. 2012). Thus, after an effective increase in 
glacier influence following major loss of ice masses, 
distinct shifts in the water supply are expected (Milner 
et al. 2009). This will cause temporary (particularly in 
summer) but also permanent desiccation within riverine 
sectors due to the resulting absence of direct (melt water 
effluent) and indirect (replenishing groundwater storage) 
glacier influences (Tockner and Malard 2003). 

In glacier-fed rivers, the relatively long history of 
research targeting the macrozoobenthos has produced 
important insights, such that its community patterns 
(e.g., species composition and abundances, and ecologi- 
cal traits) have been related to riverine habitat condi- 
tions subject to the varying influences of the glaciers 
that supply them (Ilg and Castella 2006; Fiireder 2007; 
Brown et al. 2010; Finn et al. 2010; Muhlfeld et al. 2011; 
Jacobsen and Dangles 2012). Thus, changes in climate 
and glacierization are expected to result in correspond- 
ing changes, including species extinction, in the macro- 
zoobenthos of distinctly glacier-driven habitats. For 
other organismal groups comprising the benthic biota of 
these habitats, our limited knowledge has hindered an 
understanding of their relationships with abiotic and 
biotic characters and of the expected alterations thereof 
in response to ongoing climate change (Milner et al. 
2009). This deficit derives from the widely substrate- 
biased observations of algae (epilithic; Rott et al. 2006) 
and fungi (litter traps; Gessner and Robinson 2003) and 
the paucity of studies of other groups such as bacteria 
(Battin et al. 2001, 2004) and the meiofauna (Eisendle 
2008). 



Among the meiobenthos, nematodes act as a relatively 
diverse colonizer at the glacier source and contribute, 
both in terms of abundance and diversity, to the distinct 
habitat templates further downstream (e.g., proglacial, 
glacio-rhithral) (Eisendle 2008). These findings suggest 
the important participation of nematodes as mediators in 
nutrient cycling processes in glacier-fed rivers, as similar 
to that presumed for algae as primary producers and for 
bacteria and fungi as decomposers. This role as essential 
nutrient mediators is a function of both the broad 
resource usage of nematodes and their many predators 
(Traunspurger 2002; Beier et al. 2004; Muschiol et al. 
2008; Spieth et al. 2011). However, it can also be assumed 
that nematodes, and members of the meiobenthos in gen- 
eral, are contributing factors to the as yet poorly under- 
stood resources of macrozoobenthic consumers in glacier 
rivers (Fiireder et al. 2003). Consequently, the involve- 
ment of the meiobenthos in food-web relationships and 
ecosystem performance must be considered in any discus- 
sion of the changes predicted for glacier-fed rivers and 
their biota. 

In an unprecedented approach, this study investigated 
the hitherto widely neglected micro- and meiobenthos 
(MMB: bacteria, fungi, algae, protists, nematodes, rotifers, 
copepods, tardigrades) in the Moll River catchment (MC) 
and the Kleinelendbach stream catchment (KC), both 
located in the Austrian Alps. These two differently glaci- 
ated catchments are representatives of past and ongoing 
glacier retreat patterns. Thus, in the MC, the large abla- 
tion zone of the dominant Pasterze glacier has been rap- 
idly retreating whereas the slow retreat of the KC glaciers 
over the last few decades reflects previous major losses of 
their ablation zone. Overall, the MC environment is con- 
sidered to be harsher than that of the KC mainly due to 
the extent of glacierization and the rapidly retreating ice 
masses. The two catchments have established riverine 
reaches that greatly differ in terms of harshness, a conse- 
quence of the spatiotemporally varying influences of their 
respective glaciers. Accordingly, we hypothesized: (1) 
significant differences in the MMB between the 
catchments, their reaches, and the relevant seasons 
(before, during, and after the glacier melt), and (2) that 
the changing habitat conditions produced by the varying 
glacier influences would be accompanied by significant 
alterations in the MMB. 

Material and Methods 
Study sites 

The MC and the KC (Nationalpark Hohe Tauern, Eastern 
Alps, Austria; Figs. 1, 2) were investigated at the beginning 
(GMI), middle (GMII), and end (GMIII) of the glacier 
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Figure 1. Photographs of the Moll River (MC) and Kleinelendbach 
stream (KC) catchments. 

melting period in 2010. The two catchments mainly differ 
in their areas, degree of glacierization, number of glaciers, 
and glacial retreat patterns, both during the past several 
decades and throughout the study period (Table 1). These 
differences imply different stages of succession and the 
distinct habitat conditions of the investigated river reaches 
(Table 2), which cover the youngest as well as the oldest 
potentially accessible sectors with respect to present-day 
glacial retreat and height during the Little Ice Age (1850 
AD). At each catchment, a water gauge was installed to 
measure water temperature (°C), pH, electrical conductiv- 
ity (fiS cm -1 ), and turbidity (NTU) (MS5, Otthydromet, 
Kempten, Germany). Discharge data for the MC were 
obtained from the Hydropower Agency and water level 
data at the KC from an additional gauge (Data logger, 
Orpheus Mini, Otthydromet). The vegetation of the two 
catchments consists mainly of alpine meadows. 

Field work 

Habitat stability was evaluated by the bottom component 
of the Pfankuch index (BI [bottom stability], Pfankuch 
1975) for each date and reach. Point measurements (water 
temperature [°C], pH, electrical conductivity [/.(S cm -1 ], 
0 2 [%]; measuring equipment: HI-98107pHep, Hanna 
Instruments and Multi 3420, WTW, Weilheim, Germany) 



and PVC corer samples (benthic microbiota [n = 4; sam- 
pling area: 4 cm 2 ], fungi [n = 5; 12.56 cm 2 ], protists and 
algae [n = 2; 12.56 cm 2 ], and invertebrates [n = 5; 
12.56 cm 2 ] of the upper river bed sediments [0-10 cm]) 
were obtained for each date and reach, except MC3 and 
MC4. At the latter, because of the homogenously silted 
habitat, a smaller corer (4 cm 2 ) was used for sampling 
each organismal group. Additionally, 6-10 stones were 
brushed for diatom samples (Rott et al. 2006), except at 
MC3 and MC4, where about 1 cm 2 of the uppermost 
epipelic substrate was removed with a spoon. After sam- 
pling, water depth (WD, m) and water velocity (WV, 
m -sec. Q]y|]j[ 3350, Fm. Greisinger, Regenstauf, Germany) 
were determined for the corer samples. Organismal groups 
were processed according to the methods listed in Table 
SI. 

Modified glaciality index (Gl m ) 

As in previous studies of glacial river macroinvertebrates 
by other authors (Ilg and Castella 2006; Brown et al. 
2010), we calculated a glaciality index (GI), albeit in a 
modified version (GI m ). The GI m integrates three habitat 
characteristics that change with changing glacier influence, 
namely, water temperature (T), electrical conductivity 
(EC), and BI. The reciprocal of the BI (BI re2 ) was yet used 
in the GI m to ensure that all parameters therein decrease 
with increasing glacier influence. The three variables were 
standardized before running a noncentric principal com- 
ponents analysis (NPCA), whose first- axis ordination 
scores yielded the GI m values for each sampling occasion 
(each reach at each date). 

Community characters 

Microbiota morphotypes were classified into cocci, rods, 
unbranched filamentous bacteria (UFB) (size classes of 
<10 and > 10-20 /mi), cyanobacterial colonies and 
threads, and heterotrophic and autotrophic flagellates 
(size classes of <3, >3-6, >6-10, >10 /im). Fungal CFU 
(colony-forming units) counts were measured as the total 
for each replicate and for the dominant species on each 
sampling occasion. Fungal species numbers were deter- 
mined in total for each reach. Abundances were defined 
at the genus level for protists and for algae in the Bouinse 
samples. Nematode (five replicates) and rotifer (three to 
five replicates) species abundances were determined 
whenever possible; invertebrate abundances were deter- 
mined based on coarse taxonomic resolution (five repli- 
cates). 

Nematode maturity was calculated as MI = £ V; x pi 
(Bongers 1990), where V; is the c-p value, and p ; the fre- 
quency of taxon i. The Shannon diversity of nematodes 
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Figure 2. Studied river reaches (MC1-MC4, KC1-KC4) within the two catchments (MC, Moll River catchment; KC, Kleinelendbach stream 
catchment) in the Eastern Alps (Austria). 



Table 1. Characteristics of the Moll River (MC) and Kleinelendbach 
stream (KC) catchments. 





MC 


KC 


Geographic coordinates (N/E) 


47°04'/12°44' 


47°04'/13°16' 


Catchment area (km 2 ) 


36 


12 


Number of glaciers 


9 


2 


Glacier area (km 2 ) 


22 


3 


Percentage glacierization 


61 


25 


Main glacier 


Pasterze glacier 


Kleinelendkees 






glacier 


Retreat 2010 (m) 


40 


2.7 


Retreat since 1998 (m) 


403 (debris-free 


26 




tongue) 





and rotifers was calculated as H' -£} p ; log (p ; ), where p 
is the frequency of taxa;. Diatom diversity is given in 
species numbers. 

Data analysis 

We first tested whether catchment, reach, and season or a 
combination thereof accounted for significant differences 



in the community structure of each group of organisms 
(microbiota morphotype, fungal CFU counts, algae/pro- 
tists complex, invertebrates, and nematode and rotifer 
species). As the respective data did not meet parametric 
requirements (normal distribution, variance homogene- 
ity), a nonparametric three-way analysis of variance 
(PERMANOVA; factors: catchment, reach, season) was 
applied based on the Bray Curtis resemblance matrix of 
the square-root-transformed data of each organismal 
group. The null hypothesis (no significant influence of 
any factor or of any interaction between them) was 
rejected at P < 0.01 (9999 permutations). To avoid nega- 
tive estimates of variation, we pooled the factor "catch- 
ment" for certain groups (microbiota, invertebrate coarse 
level) within the PERMANOVA routine (Anderson 2001). 
The PERMANOVA settings were as follows: partial sums 
of squares, reduced model for residuals permutation 
(except nondiatom algae and protists, and rotifers: unre- 
stricted permutation due to replicate number), fixed 
effects sum to zero for mixed terms and a nested design 
with catchment and date as fixed factors, and reach factor 
nested in catchment. 
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Table 2. Abiotic characters of the river reaches: altitude (ALT, m asl.), distance to uppermost glacier margin (DG, km), reach age (RA, years since 
deglacierization; data according to Lang [1993] for KC and according to Lieb and Slupetzky [2004] for MC), and stream width (SW, m). Point 
measurement ranges during the sample period 2010 were determined for water temperature (T, °C), electrical conductivity (EC, /iS cm -1 ), pH, 
oxygen (0 2 , %), and the bottom index (Bl; according to Pfankuch [1975]). Ranges and averages of gauge measurements (in parentheses) at KC 
and MC from July to August 2010: T, EC, and pH; 0 2 is replaced by turbidity (*NTU) for the gauge measurements. WD, VW = mean water depth 
(m) and velocity (m sec -1 ); n = 22 for MC1 and 45 for other reaches. Gl m , mean values of the modified glaciality index for each reach. 





ALT 


DG 


RA 


SW 


T 


pH 




EC 


Bl 


o 2 




WD 


WV 


Gl m 


KC1 


2241 


0.5 


35 


3-6 


7.9-9.3 


7.9- 


-8.3 


6-70 


44 


93-105 




0.16 


0.09 


-0.51 


KC2 


2190 


1.9 


70 


3-10 


4.3-9.4 


7.4- 


8 


10-65 


33 


99-102 




0.16 


0.08 


-0.11 


KC3 


2112 


1.9 


-160 


3-10 


5-7.5 


8.1 1 


5-8.23 


8-38 


33 


95-104 




0.23 


0.16 


-0.27 


KC4 


2112 


1.2 


-160 


0.5-1 


8.3-11.8 


7.3- 


-7.6 


74-120 


21 


98-1 04 




0.09 


0.01 


1.24 


MC1 


2100 


0 


0 


1-3 


0.3-0.5 


8.1- 


-8.3 


5-28 


58 


97-98 




0.10 


0.23 


-1.54 


MC2 


2070 


0.9 


30 


3-7 


0.6-1.4 


7.9- 


-8.6 


28-90 


51 


97-98 




0.47 


0.07 


-0.88 


MC3 


2070 


1.2 


48 


0.5-1 


8.8-14.5 


8.1- 


-8.8 


1 1 5-267 


25 


98-99 




0.19 


0.002 


1.04 


MC4 


2070 


1.3 


54 


0.5-1 


6-15.4 


7.9- 


-8.05 


64-270 


25 


97-101 




0.13 


0.0 


1.02 


KC 










0.3-12.5 (6.3) 


6.9- 


-7.8 (7.4) 


9-34 (23) 




*NTU 1 


369 (4) 








MC 










0.03-6.1 (4.4) 


6.2- 


-11.0 (8.9) 


0-129 (34) 




*NTU 1 


3000 (515) 









The GI m fulfilled parametric requirements (normal dis- 
tribution; homogeneity of variances) and a one-way 
ANOVA (P < 0.01) was used to test for significant differ- 
ences between the reaches, with a post hoc test (Tukey- 
HSD-Test; P < 0.05) used to identify distinct groupings 
therein. A potential relationship between habitat-specific 
glacial influence (GI m ) and community parameters (abun- 
dances, taxa richness, biotic indices) was evaluated by 
regression analysis. This step was preceded by using 
correlation patterns to determine the linearity between 
raw and transformed organism data versus the GI m . The 
parametric requirements of the residuals (normality, 
homoscedasticity) were first tested by graphical evaluation, 
followed by a Breusch-Pagan test to exactly evaluate devia- 
tions from homoscedasticity. These analyses confirmed the 
validity of using untransformed taxonomic richness (all 
groups except protists) and loglO (x + 1) -transformed 
abundance data, with the exception of the reciprocals of 
rotifer mean abundances in the regression analysis. 

The evaluation of potential relations between MMB 
community structure and abiotic habitat conditions was 
complemented by correlation analyses using Spearman's 
rho to test for: (1) relationships between stream patchi- 
ness descriptors (WV, WD) and community characters 
(abundances, taxonomic richness, and biotic indices), and 
(2) relations between the resemblance matrices of the 
complex community structure of all organismal groups 
(except epilithic diatoms) and basic habitat factors (T°, 
EC, Bl, mean WT, and WV) as determined by BIOENV 
analysis. Therein, the Bray-Curtis matrix of the log 10 
(x + 1) -transformed averages of all organisms groupings 
(except epilithic diatoms) of each sampling occasion was 
used together with the abiotic matrix (Euclidean distance) 
including hydrochemistry, mean WD and WV, and the Bl 
after normalization (WV, WD, and conductivity were 



square-root- transformed prior to normalization). The 
BIOENV output yields a global significance level 
for Spearman's rho and the variables involved in its 
calculation. 

Statistics were performed with SPSS (descriptive part, 
Spearman rank correlation, Chicago, IL), PRIMER 6.1.13 
and PERMANOVA 1.0.3 (PRIMER LTD 2009, Plymoth, 
U.K.). 

Results 

Abiotic characters 

Analyses of water temperature and turbidity indicated a 
minor influence of glacial ablation for the main channel 
of the KC compared to that of the MC (Table 2). Based 
on water temperature and Bl, extremely harsh conditions 
were determined to prevail in the proglacial complex 
(MCI, MC2), harsh conditions along the KC main chan- 
nel sites (KC1-KC3), more favorable (benign) conditions 
for the sandur complex (MC3 and MC4), and the least 
harsh conditions at KC4 (Table 2). There was no obvious 
longitudinal trend of habitat amelioration along the KC 
main channel. The site-specific harshness or absence 
thereof was also well reflected by the GI m (Fig. 3C; 
Table 2). Representing the NPCA first-axis scores it 
explained 70% of the variability, with higher loadings for 
BI rez (0.939), and water temperature (0.854) than for EC 
(0.686). The GI m significantly differed between the inves- 
tigated reaches (ANOVA: F = 30.8, P < 0.01 df = 7 ), 
distinguishing, with only one overlap, three groups of 
river reaches (post hoc Tukey test; P < 0.05): the MC 
proglacial (MCI and MC2), the main channel reaches 
(MC2, KC1-KC3), and the low-flow reaches (MC3, MC4, 
and KC4). 
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(A) 



KC1 KC2 KC3 K.C4 MCI MC2 MC3 MC4 



(B) 



(C) ,, M 




Kfr K?2 KC~1 K?4 M^i hiC2 tJci M<^4 



Figure 3. Spatiotemporal distribution of (A) counts of microbiota as 
(o, bacteria; □, unbranched threads; x, cyanobacteria; A, HNF; +, ANF 
- all Counts/cm 2 /mL; 0, fungal CFU); (B) abundances of dominant 
invertebrate groups (o, nematodes; □, rotifers; x, copepods; A, 
tardigrades; lnd/10 cm 2 /ml_); and (C) the Gl m . 



Overall catchment and reach observations 

In total, the KC harbored more species than the MC, 
which was characterized by less diverse assemblages for 
all groups but diatoms (Table 3). More than double the 
number of invertebrate and algal (other than diatoms) 
taxa were present in the KC than in the MC, while the 
number of protist genera nearly tripled from the MC to 
the KC. Species numbers were highest for diatoms, 
followed by nematodes and rotifers, and lowest for 
protists. The GI m reach-specific grouping was better 
reflected by nematode maturity than by other diversity 
measures (Table 3), which were nonetheless lower in the 
main channel reaches of the MC than in those of the 
KC, where they did not consistently increase down- 
stream. Most of the diatom (MC3) and nematode 
(MC4) species inhabited the sandur complex. The lowest 
total taxonomic richness was not observed at the youn- 
gest (MCI), but at the second youngest reach (MC2), 
which although nearly the same age as KC1 had half the 
number of taxa than KC1. The difference between the 
similarly aged KC3 and KC4 was also considerable but 
not as large. 

Eumonhystera filiformis (Bastian, 1865), Mononchus 
truncatus (Bastian, 1865), and M. sandur (Eisendle 2008) 
dominated the MC floodplain sites. E. pseudobulbosa 
(Daday, 1896) and E. longicaudatula (Gerlach and Rieman 
1973) were dominant in the habitats of the main channels 
of the KC, and E. vulgaris (De Man, 1880), and two Tob- 
rilidae species in the KC4 habitat, including Tobrilus ap- 
ophysis (Steiner, 1919), also being highly abundant at 
KC3. Bdelloid rotifers dominated the rotifer community 
at all sites. Mucor sp., Fusarium sp., and Cladosporium sp. 
dominated the fungal community. Achnanthidium minu- 
tissimum (Kiitzing) Czarnecki was the dominant species 
among diatoms. A comprehensive taxon list is provided 
in Table S3. 



Microbiota 

Bacteria were dominated by cocci. Among the flagellates, 
autotrophic forms prevailed in most of the replicates 
(93%). The ranking of mean abundances among reaches 
varied among the microbiota groups (Table 4). Lowest 
means at the recently emerged MCI were determined 
only for bacteria (summed cocci and rods) and UFB. The 
mean abundances of HNF and ANF (autotrophic nanofla- 
gellates) were lowest at the main channel reaches of the 
KC and increased toward the MC proglacial 
(KCKKC3<KC2<MCKMC2). Bacteria, UFB, and flagel- 
late (ANF and HNF) abundances were high and followed 
similar patterns of MC4>MC3>KC4, unlike the abun- 
dances of cyanobacteria (MC4>MC1>KC1) and fungi 
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Table 3. Taxon numbers summed for each reach and catchment as genus numbers for algae other than diatoms (AlgO; incl. cyanobacteria) and 
protists (PRO) and species numbers for other organisms (FUN, fungi; DIA, diatoms; ROT, rotifers; NEM, nematodes). The genera numbers of the 
latter are given in parentheses. Ml, nematode maturity; INV, higher invertebrate taxa level; H', Shannon diversity for the respective groups; SN, 
sums of the taxa for each reach and catchment. 





FUN 


DIA 


AlgO 


PRO 


ROT 


H' ROT 


NEM 


FT NEM 


Ml 


INV 


SN Tot 


KC1 


14(9) 


46 (24) 


14 


5 


12 (6) 


0.1 


27 (19) 


2.4 


2.3 


12 


130 


KC2 


15(9) 


60 (29) 


15 


7 


23 (9) 


0.4 


34 (20) 


2.3 


2.3 


13 


167 


KC3 


14(7) 


35 (24) 


14 


6 


19(9) 


0.3 


39 (25) 


2.5 


2.2 


12 


139 


KC4 


15(9) 


71 (33) 


19 


10 


27 (12) 


1.7 


40 (21) 


2.7 


2.9 


13 


195 


MC1 


15 (10) 


27 (18) 


1 


2 


6(3) 


0.1 


18 (12) 


2.6 


2.0 


3 


72 


MC2 


12 (7) 


26 (14) 


1 


1 


1 (D 


0 


19(13) 


2.1 


2.0 


4 


64 


MC3 


12 (8) 


91 (36) 


6 


3 


20 (9) 


1.1 


37 (21) 


1.6 


2.6 


3 


172 


MC4 


9(6) 


84 (36) 


5 


1 


17 (8) 


1.2 


44 (23) 


1.5 


2.6 


6 


166 


KC 


29 (13) 


1 08 (39) 


20 


11 


39 (14) 


1.5 


69 (34) 


2.9 


2.8 


15 


291 


MC 


25 (15) 


142 (42) 


9 


4 


26 (10) 


1.0 


57 (30) 


1.6 


2.6 


7 


270 



(KC4>MC3>MC2) (Table 4). The seasonal patterns of 
total abundances also differed between these groups 
although there was no consistent decrease among and 
between them at either reach during the summer ablation 
period GMII (Fig. 3A). The spatiotemporality of the bac- 
teria and nanoflagellate complex was significantly influ- 
enced by reach and by the interaction between reach and 
time, but not by catchment, which had a slightly greater 
influence than either of the other two factors on the vari- 
ability in the algae/protist complex (Table 5). 

Invertebrates 

Rotifers and nematodes dominated the invertebrate 
fauna. Copepods and macroinvertebrates such as chir- 
onomids or planarians (data not shown) were extremely 
scarce not only at all MC reaches but also along the 
main channel of the KC (Table 4). Copepods were only 
present with considerable abundances at KC4. Nema- 
tode and rotifer abundances were lowest at the harsh 
proglacial (MC1<MC2), whereas copepod abundances 
were lowest at MCI and KC3 and tardigrades ranked 
totally differently (KC2<KC3<KC1). The rankings of the 
highest abundance differed between groups and could 
not be consistently ascribed to either deglacierization 
stage or habitat conditions: MC4>KC4>MC3 (nema- 
todes), KC4>MC4>MC3 (rotifers), and KC4>KC1>KC2 
(copepods). Spatiotemporal abundances did not undergo 
a uniform decline in summer (Fig. 3B). Analysis of the 
complex community structures by resemblance matrices 
within the PERMANOVA analysis showed that while 
for nematodes reach and an interaction between reach 
and date were relevant, only reach explained the signifi- 
cant differences in rotifer resemblance patterns; catch- 
ment had no significant influence on either group 
(Table 5). 



Relation between abiotics and community 
patterns 

The relations between the GI m and the community char- 
acters that fulfilled the regression requirements (18) are 
listed in Table 6. The GI m explained >70% of the varia- 
tions in nematode maturity and abundances and >50% of 
the variations in six other characters (e.g., nematode, roti- 
fer, and diatom taxonomic richness). Six factors did not 
significantly alter in parallel with changes in the GI m . WD 
correlated negatively with tardigrades, HNF (both 
P < 0.01), nematodes, and copepods (both P < 0.05); 
WV negatively with nematodes, tardigrades, bacteria, 
threads, HNF, ANF (all: P < 0.01), rotifers, and cyano- 
bacteria (both P < 0.05). A significant positive correlation 
was established among invertebrates groups and among 
microbiota morphotypes (P < 0.01) (all Spearman's rho). 
Finally, according to the BIOENV, the variability in BI 
alone was most strongly related to the complex variability 
in the MMB resemblance patterns (Spearman's rho: 0.54; 
P = 0.01, number of permuted rho's >global rho = 0). 

Discussion 
General aspects 

The distinctly different past and present characteristics 
and activities of the studied glaciers at the two catch- 
ments manifested as distinct habitat templates within a 
range of extremely rough to benign conditions. These 
conditions were condensed into three groups by the GI m : 
extremely harsh (MCI and MC2), harsh (KC1-KC3), and 
benign (MC3, MC4, and KC4). The MC is representative 
of a high degree of glacierization, but also of a glacier 
(Pasterze glacier) that has rapidly decreased in size during 
the last few decades. The latter accounts for the two 
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distinct sections of the MC: (1) a highly unstable progla- 
cial, comprising the glacier source (MCI) and MC2, 
which consists mainly of unconsolidated moraine mate- 
rial, and (2) the consecutive silted floodplain area (Krain- 
er and Poscher 1992). The ongoing rapid loss of ice 
masses during the study period (Fischer 2010) led to a 
strong glacier influence at the proglacial (e.g., low water 
temperature and a low GI m ) that was distinctly modified 
along the silted floodplain area (e.g., high water tempera- 
ture and high GI m ). Such silted floodplain areas can often 
be ascribed to large melt water streams (Krigstom 1962; 
Hodgkins et al. 2009). 

In contrast, the KC is a typical example of a less glaci- 
erized catchment where the accumulation area now pre- 
dominates after the wide and effective retreat of the 
ablation area. It is also representative of the fragmentation 
of a once contiguous glacial surface, as is similarly 
expected for many glaciers in response to predicted 
climate changes (Jiskoot and Miiller 2012). This fragmen- 
tation, resulting partly in cliff glacier remnants, was 
responsible for habitat conditions that contradicted a 
longitudinal gradient of amelioration along the KC main 
channel sites, particularly during the melting period. Cliff 
remnants at KC3 intruded upon a proposed amelioration 
gradient, as their summer-time effluents caused high cur- 
rent-water velocities at KC3. These findings suggest that 
the downstream development of glacier streams as a con- 
sequence of ongoing glacier retreat will abrogate any 
downstream longitudinal increase in benignity in glacier 
rivers (see Milner et al. 2001) when a stream accompany- 
ing the glacier remnants develops in parallel. 

One of the basic assumptions in this study, namely, 
that the two catchments cause significant differences in 
the MMB because of differences in their harshness - 
defined according to the percentage of glacierization, gla- 
cier retreat patterns, and stream ages - was not 
confirmed. Nevertheless, single community descriptors 
(e.g., algal, protist, nematode, and rotifer species num- 
bers; rotifer and nematode H' diversity; and nematode 
maturity), might better capture the harshness of the 
catchments, as evidenced by the lower values in the MC, 
indicative of higher overall catchment harshness. Struc- 
tural alterations in the macrozoobenthos were shown to 
occur with decreasing glacerization, with an increase in 
the numbers of individuals and species with ongoing 
glacier retreat in alpine areas (Fiireder 2007). Our study 
comparing the KC and MC is the first to analyse similar 
causalities for species number of the micro- and meiozoo- 
benthos, comprising protists, rotifers, and nematodes. 
Thus, as for the macrofauna, we also found evidence of 
an increase in the taxonomic richness with decreasing 
glacierization, at least to the extent occurring between the 
two catchments (from 60% to 25%). 
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Table 5. Output of the PERMANOVA for the respective groups: Pseudo-F values are shown for the factors reach (R), date (D), and catchment 
(CA) and their interactions. Values in italics indicate the contribution of estimated variations derived from the respective factor. 





R 


R x D 


D 


CA 


CA x D 


CA 


CA x D 


R 


R x D 


D 


Res 


Microbiota 


27 9** 


2.4** 


1.7 


Pooled 


1.2 




3 


32 


12 


5 


21 


Fungi CFU 


7.0** 


3 2** 


1.3 


1.1 


4.5 


2 


16 


13 


14 


3 


21 


Algae/protists 


2 3** 


1.2 


2 7** 


3.3** 


1.2 


26 


17 


27 


15 


17 


44 


Invertebrates 


33.6** 


2.0** 


2.3 


Pooled 


1.4 




5 


35 


11 


6 


24 


Nematodes 


1 2 9** 


2"]** 


1.6 


1.7 


1.4 


17 


9 


39 


21 


8 


42 


Rotifers 


11.6** 


1.5 


2.4 


1.5 


1.5 


11 


7 


30 


11 


8 


30 



Significance in the rejection of HO at **P < 0.01. Number of permutations = 9999. 



Table 6. Results of the regression analysis (R 2 , F, bO, and fo1 as the 
regression coefficients of the linear functions). 





R 2 


F 


bO; 61 




Nematode Ml 


0.74** 


63.78 


2.36 


0.29 


Nematode TR 


0.53** 


23.40 


17.18 


4.72 


Nematode H' 


2.E-03 0 


0.04 


1.88 


-0.02 


Nematode M_L10 


0.77** 


71.09 


1.94 


0.84 


Rotifera TR 


0.62** 


34.00 


9.09 


4.71 


Rotifera H' 


0.67** 


41.99 


0.59 


0.49 


Rotifera M_REZ 


0.13° 


3.32 


0.01 


-0.01 


Bacteria M_L10 


0.53** 


23.23 


6.68 


0.54 


HNF M_L10 


0.42** 


14.96 


3.01 


0.53 


Diatoms TR 


0.64** 


36.69 


30.50 


13.69 


Algae TR 


0.27* 


7.77 


6.78 


2.71 


ANF M_L10 


0.23* 


6.41 


3.76 


0.41 


Cyanobacteria M_L10 


8.E-08 0 


0.00 


2.48 


-0.02 


UFB M_L10 


0.62** 


34.25 


4.09 


0.69 


Algae M_L10 


0.29** 


8.58 


4.20 


0.95 


Protists M_L1 0 


0.02° 


0.46 


1.63 


0.13 


Protists TR M_L10 


0.11° 


2.69 


0.75 


0.45 


Fungal CFU M_L10 


0.12° 


2.97 


2.19 


0.15 



Ml, maturity; TR, taxonomic richness; H', Shannon diversity; M, mean 
abundance; L10, log 10 -transformed data; REZ, reciprocal value; UFB, 
unbranched filamentous bacteria. 

**Highly significant: P<0.01, 'significant: P < 0.05 and °no signifi- 
cant relation between the Gl m and the respective community parame- 
ters. 

Between-catchment comparisons, as discussed above, 
are generally rare in glacier river studies. More often, the 
focus is on the distribution of community patterns (abun- 
dances, richness) as a function of varying habitat harsh- 
ness along three main axes: (1) the longitudinal decrease 
in harshness with increasing distance from the glacier 
margin, which should typically be paralleled by an 
increase in stream age; (2) a seasonal decrease in harsh- 
ness outside the glacier ablation period (1 and 2: Ward 
1994; Milner et al. 2001); and (3) a lateral component 
enabling the existence of less extreme rhithral, glacio- 
rhithral, and krenal habitats (Tockner and Malard 2003; 
Uehlinger et al. 2003). Most of the studies referring to 



these axes have reported an increase in structural com- 
plexity with respect to the diversity and/or abundances of 
their studied organisms together with a decreases in 
harshness (Battin et al. 2001; Lods-Crozet et al. 2001a,b; 
Robinson et al. 2001; Burgherr et al. 2003; Biirgi et al. 
2003; Malard et al. 2006; Eisendle 2008). However, obser- 
vations refuting these basic schemes have also been made, 
for example, a report of a decrease in fungal species num- 
bers in the river downstream (Gessner and Robinson 
2003) and the results of this study (but see below). 

In this study, taxon numbers of benthic fungi and of 
most of the other groups decreased (see Table 3) down- 
stream from the MC glacier source to MC2, and from 
KC2 to KC3 although by different amounts. MMB abun- 
dances also partly decreased in the river downstream 
along the MC and KC main channel sites. These observa- 
tions deviate from the abovementioned scheme in glacier 
river research, in which a longitudinal downstream 
increase in taxonomic richness and abundances in rela- 
tion to downstream habitat amelioration and succession 
progress (stream age). However, these deviations are 
consistent with the particular catchment characters in this 
study, as already discussed herein, namely, the instability 
of MC2 and the lengthened, glacier-accompanied main 
channel course of the KC. Finally, to answer the question 
whether harshness or stream age has greater determining 
relevance for the MMB, comparisons among KC3 and 
KC4, as well as KC1 and MC2, and KC4 and the sandur 
complex lead to the conclusion that benignity and harsh- 
ness are more important than stream age, both within 
and between catchments. 

In its original form, the GI was introduced to improve 
the relational analysis of glacier influence on macrozoo- 
benthic patterns (Ilg and Castella 2006). The GI integrates 
within-stream conditions that change with varying glacier 
influence (water temperature, EC, BI, and turbidity). 
Despite several shortcomings, the index is relatively easy 
to apply and, in most cases it distinguishes complex habi- 
tat conditions, relating them, more or less strongly, to 
community patterns (Ilg and Castella 2006; Brown et al. 
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2010). The reasonable reach grouping confirmed the util- 
ity of the modified version developed for this study. More 
importantly, perhaps, ours is the first application of the 
GI to the MMB and it revealed well-established relation- 
ships for most of the micro- and meiobiota characters. 
These relationships were often stronger than those previ- 
ously reported for the few macroinvertebrate structural 
parameters (richness, abundances) examined thus far 
(Brown et al. 2010). Specifically, the GI m clearly indicated 
the positive response of certain community parameters 
(e.g., taxonomic richness of nematodes, diatoms, and 
other algae; abundances of bacteria, flagellates, algae, and 
nematodes) to decreasing glacier influence, which in itself 
entails complex and reach specifically altered conditions 
between and among catchments. 

Predictions on the nature of the MMB changes along 
glacier-dependent river systems, however, should be gen- 
eralized only with caution, based on the varying results 
described herein for the different groups. A summary of 
the results obtained for relation patterns between the 
MMB and several abiotic parameters clearly showed, for 
example, that an increase in water temperature is insuffi- 
cient to promote an increase in the abundance or taxo- 
nomic richness of most groups when it is accompanied 
by high-current velocity and low-bed instability. The pre- 
dicted summer-time shift to precipitation-driven effluents 
in alpine catchments as discussed in Milner et al. (2009) 
might therefore have similar repressive effects on MMB 
richness and density, due to their flashy, erosive, and less 
predictable run-off character. A shift of water sources 
might also bear serious consequences for catchment func- 
tionality if there is a loss of benign habitats in response 
to a loss of groundwater recharge function during the 
ablation period (Tockner and Malard 2003). It is particu- 
larly these habitats that, regardless of their bed sediment 
texture, are important local hot spots of diversity and 
thus of functionality. 

Selected community aspects 

This study adds new and important information regard- 
ing the autotrophic component of glacier-fed rivers. Our 
results suggest that high-diatom species richness is not 
restricted to the epilithon or to those habitats with 
diminished glacial influence and are thus in contrast to 
previous studies (see Rott et al. 2006). In particular, the 
number of diatom species at silted floodplains highlights 
the need for greater attention to these habitats in future 
research. However, it is not only diatom species numbers 
that seem to be distinctly underrepresented when exclu- 
sively determined on epilithic substrates, as indicated by 
other algal species number among KC reaches (compare 
Biirgi et al. 2003). A comparison of the photoautotrophic 



components among the MC and KC reaches showed the 
delayed succession of a diverse autotrophic community 
not limited to diatoms in riverine habitats shaped by 
rapid deglacierization. Finally, ANFs, although widely 
neglected, may be important primary producers in 
glacier-fed river habitats. As a subject of cryoconite 
research, their general ecological importance within gla- 
cier biomes has already been demonstrated (Hodson et al. 
2008; Anesio et al. 2009). Based on their distribution pat- 
terns, as determined in our study, ANFs can be consid- 
ered as a constant and seemingly adaptive component 
within carbon cycling processes in glacier-fed river habi- 
tats driven by varying glacier influences. 

Despite consistently strong differences in the thermal 
and hydrological regimes of the two catchments, there 
were no significant corresponding differences in the 
resemblance patterns of the bacterial and flagellate micro- 
biota. Specifically, microbiota densities might be positively 
rather than negatively affected by a huge loss of ice 
masses. This conclusion is based on the abovementioned 
growing appreciation of glaciers as important habitats of 
microbial organisms, thereby fulfilling important ecologi- 
cal functions on a regional as well as a global scale 
(Kastovska et al. 2007; Mindl et al. 2007; Hodson et al. 
2008; Anesio et al. 2009; Stibal et al. 2012; Viles 2012). 
Bacteria have been shown to rapidly alter their genetic 
expression patterns, changing from swimming to sessile 
life forms (Hall-Stoodley et al. 2004). This could be 
advantageous in shifts from glacial to riverine benthic 
habitats, such as by facilitating colonization at recently 
deglacierized sites, in favoring downstream drift, and in 
the recolonization of unstable habitats. 

The distribution patterns of the fungal CFU counts and 
fungal species numbers did not show a uniformly increas- 
ing trend in the shift from harsh to benign conditions or 
from high to low glacierization, suggesting their indepen- 
dence from varying glacial influence and thus from 
glacially driven hydrological and thermal regimes. The taxa 
from the benthic substrate clearly differed from those 
described in glacier river studies on fungi, in which leaf- 
litter traps were used (Gessner and Robinson 2003). The 
MC and KC taxa partly resembled the spectra of terrestrial 
taxa related to understory vegetation (De Bellis et al. 
2007), corroborating the assumed role of fungi in connect- 
ing aquatic habitats and their catchments (Jobard et al. 

2010) but also implying changes in the benthic fungal 
community in response to those in the surrounding vege- 
tation (Walther et al. 2005; Parmesan 2006; Leonelli et al. 

2011) . The observed fungal patterns might reflect the 
importance of fungi as an essential nutrient resource in 
high-altitude reaches, particularly for smaller invertebrates 
(e.g., nematodes and rotifers), by providing suitable 
high-energy and high-quality "mouthfuls" of food. This, 
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together with the reliance of the fungal taxa on their sur- 
roundings, would add a further, previously unrecognized 
component to both the intrinsic relationships between 
potential resources and their consumers and the alterations 
thereof by ongoing changing environmental conditions. 

Among those invertebrates for which species composi- 
tion was investigated, habitat conditions and deglacieriza- 
tion stages contributed less to the similarity patterns of 
rotifers than of nematodes. This finding can be explained 
by the clear predominance of bdelloid rotifers, a group 
particularly well adapted to hardy conditions and with a 
number of pioneer species (Wallace et al. 2006; Hodson 
et al. 2008; Webster-Brown et al. 2010). Both the behav- 
ioral adaptations of monogononts searching for shelter in 
deeper interstitial layers (Schmid-Araya 2000) and the 
diurnal migrations of monogononts and nematodes aimed 
at high-flow avoidance (Smith and Brown 2006) would 
account for the maintenance of low, but permanent densi- 
ties of these groups despite frequent high-flow conditions. 
These activities would also provide a behavioral advantage 
for these organisms in the future in response to climate 
change, when they are potentially confronted with an 
increased glacial melt and the prolongation thereof as well 
as less predictable precipitation flood events. 

Although less abundant at particularly harsh sites, nem- 
atodes were more diverse than rotifers at each reach and 
catchment, with the highest abundances at the extremely 
differently structured benign sites. Therefore, nematodes 
might be able to well cope with a variety of reach condi- 
tions, a conclusion supported by the variety of running 
waters harboring nematodes (Eisendle 2009 and references 
therein). The nematode resemblance patterns clearly 
differentiated between the different KC and MC habitat 
types, and the nematode maturity was related to the GI m - 
established reach grouping. Furthermore, the observed 
nematode feeding types (data not shown; Yeates et al. 
1993), were in good accordance with the respective food 
resources of each investigated reach. Together, these 
results support the usefulness of nematodes as indicators 
of glacially induced habitat templates as well as the devel- 
opment of benthic resource variety along glacier-fed 
rivers. Finally, we tentatively predict the loss of Hofma- 
enneria hazenensis (Mulvey 1969) (first recorded for the 
Arctic), with the loss of glacial influence, because in this 
study it inhabited only those river sites with distinct sur- 
face glacier influences (all except KC4). 

Acknowledgments 

This study was financed entirely (incl. cooperations) by a 
research grant (Austrian Science Fund - FWF P 21922) to 
the first author. We are grateful to Cornelia and Michael 
Lunder, Eva Wilhelm, Johann Flbckner and Ambros 



Unterkircher for aid with sampling, and to M. Lunder 
and A. Unterkircher for their assistance in installing the 
water gauges. We thank the Verbund Tourismus and the 
GroglocknerhochalpenstraEen AG (GROHAG) for 
enabling free access to the toll roads. We also thank 
Erwin Hasslacher (Nationalpark Hohe Tauern) for his 
ongoing support. We thank Doris Gesierich for diatom 
determination and Lee Brown for comments on the glaci- 
ality index. 

Conflict of Interest 

None declared. 
References 

Anderson, M. ). 2001. Permutation tests for univariate or 

multivariate analysis of variance and regression. Can. J. Fish. 

Aquat. Sci. 58:626-639. 
Anesio, A. M., A. ). Hodson, A. Fritz, R. Psenner, and B. 

Sattler. 2009. High microbial activity on glaciers: importance 

to the global carbon cycle. Glob. Change Biol. 15:955-960. 
Battin, T. J., A. Wille, B. Sattler, and R. Psenner. 2001. 

Phylogenetic and functional heterogeneity of sediment 

biofilms along environmental gradients in a glacial stream. 

Appl. Environ. Microbiol. 67:799-807. 
Battin, T. J., A. Wille, R. Psenner, and A. Richter. 2004. 

Large-scale environmental controls on microbial biofilms in 

high-alpine streams. Biogeosciences 1:159-171. 
Beier, S., M. Bolley, and W. Traunspurger. 2004. Predator-prey 

interactions between Dugesia gonocephala and free-living 

nematodes. Freshw. Biol. 49:77-86. 
Bongers, T. 1990. The maturity index: an ecological measure 

of environmental disturbance based on nematode species 

composition. Oecologia 83:14-19. 
Brown, L. E., D. M. Hannah, and A. M. Milner. 2009. ARISE: 

a classification tool for Alpine River and Stream Ecosystems. 

Freshw. Biol. 54:1357-1369. 
Brown, L. E., A. M. Milner, and D. M. Hannah. 2010. 

Predicting river ecosystem response to glacial meltwater 

dynamics: a case study of quantitative water sourcing and 

glaciality index approaches. Aquat. Sci. 72:325-334. 
Burgherr, P., B. Klein, C. T. Robinson, and K. Tockner. 2003. 

Surface zoobenthos. Pp. 153-173 in ]■ V. Ward and U. 

Uehlinger, eds. Ecology of a glacial flood plain. Kluwer 

Academic Publishers, Dordrecht. 
Biirgi, H. R., P. Burgherr, and U. Uehlinger. 2003. Aquatic 

flora. Pp. 139-151 in J. V. Ward and U. Uehlinger, eds. 

Ecology of a glacial flood plain. Kluwer Academic 

Publishers, Dordrecht. 
De Bellis, T., G. Kernaghan, and P. Widden. 2007. Plant 

community influences on soil microfungal assemblages in 

boreal mixed-wood forests. Mycologia 99:356-367. 



2842 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



U. Eisendle-Flockner eta/. 



Small Benthic Biota in Glacier Catchments 



Dessler, A., and E. A. Parsons. 2010. The science and politics 

of global climate change: a guide to the debate. Cambridge 

Univ. Press, Cambridge. 211 pp. 
Dyurgerov, M., M. Meier, and D. Bahr. 2009. A new index of 

glacier area change: a tool for glacier monitoring. J. Glaciol. 

55:710-716. 

Eisendle, U. 2008. Spatiotemporal distribution of free-living 
nematodes in glacial-fed stream reaches (Hohe Tauern, 
Eastern Alps, Austria). Arct. Antarct. Alp. Res. 40:470-480. 

Eisendle, U. 2009. Distribution patterns of free-living 

nematodes in differently degraded river reaches of the River 
Wien, Austria. Nematology 11:509-521. 

Farinotti, D., S. Usselmann, M. Huss, A. Bauder, and M. 
Funk. 2012. Runoff evolution in the Swiss Alps: projections 
for selected high-alpine catchments based on ENSEMBLES 
scenarios. Hydrol. Process. 26:1909-1924. 

Finn, D. S., K. Raesaenen, and C. T. Robinson. 2010. Physical 
and biological changes to a lengthening stream gradient 
following a decade of rapid glacial recession. Glob. Change 
Biol. 16:3314-3326. 

Fischer, A. 2010. Gletscherbericht (2009/2010). Bergauf 
135:24-30. 

Fiireder, L. 2007. Life at the edge: habitat condition and 

bottom fauna of alpine running waters. Int. Rev. Hydrobiol. 
92:491-513. 

Fiireder, L., C. Welter, and J. K. Jackson. 2003. Dietary and 

stable isotope (<513 C, <5 15 N) analyses in alpine stream 

insects. Int. Rev. Hydrobiol. 88:314-331. 
Gessner, M. C, and C. T. Robinson. 2003. Aquatic 

hyphomycetes in alpine streams. Pp. 123-137 in J. V. Ward 

and U. Uehlinger, eds. Ecology of a glacial flood plain. 

Kluwer Academic Publishers, Dordrecht. 
Hall-Stoodley, L., J. W. Costerton, and P. Stoodley. 2004. 

Bacterial biofilms: from the natural environment to 

infectious diseases. Nat. Rev. Microbiol. 2:95-108. 
Hodgkins, R., R. Cooper, J. Wadham, and M. Tranter. 2009. 

The hydrology of the proglacial zone of a high-Arctic glacier 

(Finsterwalderbreen, Svalbard): atmospheric and surface 

water fluxes. J. Hydrol. 378:150-160. 
Hodson, A., A. M. Anesio, M. Tranter, A. Fountain, 

M. Osborn, and J. Priscu, et al. 2008. Glacial ecosystems. 

Ecol. Monogr. 78:41-67. 
Hook, R. 2005. Glacier melt: a review of processes and their 

modelling. Prog. Phys. Geogr. 29:362-391. 
Ilg, C, and E. Castella. 2006. Patterns of macroinvertebrate 

traits along three glacial stream continuums. Freshw. Biol. 

51:840-853. 

Jacobsen, D., and O. Dangles. 2012. Environmental harshness 
and global richness patterns in glacier-fed streams. Glob. 
Ecol. Biogeogr. 21:647-656. 

Jiskoot, H, and M. S. Miiller. 2012. Glacier fragmentation 
effects on surface energy balance and runoff: field 
measurements and distributed modelling. Hydrol. Process. 
26:1861-1875. 



Jobard, M., S. Rasconi, and T. Sime-Ngando. 2010. Diversity 
and functions of microscopic fungi: a missing component in 
pelagic food webs. Aquat. Sci. 72:255-268. 

Kastovska, K., M. Stibal, M. Sabacka, B. Cerna, H. 

Santruckova, and J. Elster. 2007. Microbial community 
structure and ecology of subglacial sediments in two 
polythermal Svalbard glaciers characterized by 
epifluorescence microscopy and PLFA. Polar Biol. 
30:277-287. 

Krainer, K., and G. Poscher. 1992. Sedimentologische 
Beobachtungen im Gletschervorfeld der Pasterze 
(Glocknergruppe, Hohe Tauern). Carinthia II 182:317-343. 

Krigstom, A. 1962. Geomorphological studies of sandur plains 
and their braided rivers in Iceland. Geogr. Ann. 44:328- 
346. 

Lang, H. 1993. Die Gletscher Karntens. Naturwissenschaftlicher 

Verein Karntens, Klagenfurt. 
Lecerf, A., and J. S. Richardson. 2010. Biodiversity-ecosystem 

function research: insights gained from streams. River Res. 

Appl. 26:45-54. 
Lieb, G. K., and H. Slupetzky. 2004. Gletscherweg Pasterze. 

Naturkundlicher Fiihrer zum Nationalpark Hohe Tauern 2 

(2. vollig neu bearbeitete Auflage), Innsbruck. 122 pp. 
Leonelli, G., M. Pelfini, U. Morra Di Cella, and V. Garavaglia. 

2011. Climate warming and the recent treeline shift in the 

European alps: the role of geomorphological factors in 

high-altitude sites. Ambio 40:264-273. 
Lods-Crozet, B., E. Castella, D. Cambin, C. Ilg, S. Knispel, and 

H. Mayor-Simeant. 2001a. Macroinvertebrate community 

structure in relation to environmental variables in a Swiss 

glacial stream. Freshw. Biol. 46:1641-1661. 
Lods-Crozet, B., V. Lencioni, J. S. Olafsson, D. L. Snook, 

G. Velle, and J. E. Brittain, et al. 2001b. Chironomid 

(Diptera: Chironomidae) communities in six European 

glacier- fed streams. Freshw. Biol. 46:1791-1809. 
Malard, F., U. Uehlinger, R. Zah, and K. Tockner. 2006. 

Flood-pulse and riverscape dynamics in a braided glacial 

river. Ecology 87:704-716. 
McGregor, G., G. E. Petts, A. M. Gurnell, and A. M. Milner. 

1995. Sensitivy of alpine stream ecosystems to climate 

change and human impacts. Aquat. Conserv. 5:233- 

247. 

Milner, A. M., J. E. Brittain, E. Castella, and G. E. Petts. 2001. 
Trends of macroinvertebrate community structure in 
glacier-fed rivers in relation to environmental conditions: a 
synthesis. Freshw. Biol. 46:1833-1847. 

Milner, A. M., L. E. Brown, and D. M. Hannah. 2009. 
Hydroecological response of river systems to shrinking 
glaciers. Hydrol. Process. 23:62-77. 

Mindl, B., A. M. Anesio, K. Meirer, A. J. Hodson, J. 
Laybourn-Parry, R. Sommaruga, et al. 2007. Factors 
influencing bacterial dynamics along a transect from 
supraglacial runoff to proglacial lakes of a high Arctic 
glacieri. FEMS Microbiol. Ecol. 59:307-317. 



2013 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



2843 



Small Benthic Biota in Glacier Catchments 



U. Eisendle-Flockner ef al. 



Muhlfeld, C. C, J. J. Giersch, F. R. Hauer, G. T. Pederson, G. 
Luikart, and D. P. Peterson, et al. 2011. Climate change 
links fate of glaciers and an endemic alpine invertebrate. 
Clim. Change 106:337-345. 

Mulvey, R. H. 1969. Soil-inhabiting nematodes of the orders 
Araeolaimida, Chromadorida, Enoplida and Monhysterida 
from the Canadian high Arctic. Can. J. Zool. 47:365-382. 

Muschiol, D., M. Markovic, I. Threis, and W. Traunpurger. 
2008. Predator-prey relationship between the cyclopoid 
copepod Diacyclops bicuspidatus and a free-living 
bacterivorous nematode. Nematology 10:55-62. 

Parker, S. M., and A. D. Huryn. 2011. Effects of natural 
disturbance on stream communities: a habitat template 
analysis of arctic headwater streams. Freshw. Biol. 
56:1342-1357. 

Parmesan, C. 2006. Ecological and evolutionary responses to 
recent climate change. Annu. Rev. Ecol. Syst. 37:637-669. 

Pelto, M., and C. Brown. 2012. Mass balance loss of Mount 
Baker, Washington glaciers 1990-2010. Hydrol. Process. 
26:2601-2607. 

Pfankuch, D. J. 1975. Stream reach inventory and channel 
stability evaluation. US Department of Agriculture Forest 
Service, Region 1, Missoula, MT. 26 pp. 

Robinson, C. T., U. Uehlinger, and M. Hieber. 2001. 

Spatio-temporal variation in macroinvertebrate assemblages 
of glacial streams in the Swiss Alps. Freshw. Biol. 
46:1663-1672. 

Rott, E., M. Cantonati, L. Fureder, and P. Pfister. 2006. 
Benthic algae in high altitude streams of the Alps - a 
neglected component of the aquatic biota. Hydrobiologia 
562:195-216. 

Schmid-Araya, J. M. 2000. Invertebrate recolonization patterns 
in the hyporheic zone of a gravel stream. Limnol. Oceanogr. 
45:1000-1005. 

Smith, F., and A. V. Brown. 2006. Effects of flow on 

meiofauna colonization in artificial streams and reference 
sites within the Illinois River, Arkansas. Hydrobiologia 
571:169-180. 

Spieth, H. R., T. Moller, C. Ptatschek, A. Kazemi-Dinan, and 
W. Traunspurger. 2011. Meiobenthos provides a food 
resource for young cyprinids. J. Fish Biol. 78:138-149. 

Stendera, S., R. Adrian, N. Bonada, M. Canedo-Argiiells, 
B. Hugueny, and K. Januschke, et al. 2012. Drivers and 
stressors of freshwater biodiversity patterns across different 
ecosystems and scales: a review. Hydrobiologia 696:1-28. 

Stibal, M., J. Telling, J. Cook, K. M. Mak, A. Hodson, and A. 
M. Anesio. 2012. Environmental controls on microbial 
abundance and activity on the Greenland ice sheet: a 
multivariate analysis approach. Microb. Ecol. 63:74-84. 



Tockner, K., and F. Malard. 2003. Channel typology. Pp. 

57-73 in J. V. Ward and U. Uehlinger, eds. Ecology of a 

glacial floodplain. Kluwer Academic Publishers, Dordrecht. 
Traunspurger, W. 2002. Nematoda. Pp. 63-104 in S. D. Rundle, 

A. L. Robertson and J. M. Schmid-Araya, eds. Freshwater 

meiofauna: biology and ecology. Backhuys Publishers, Leiden. 
Uehlinger, U., F. Malard, and J. V. Ward. 2003. Thermal 

patterns in the surface waters of a glacial river corridor (Val 

Roseg, Switzerland). Freshw. Biol. 48:284-300. 
Vanham, D. 2012. The Alps under climate change: 

implications for water management in Europe. J. Water 

Clim. Change 3:197-206. 
Viles, H. A. 2012. Microbial geomorphology: a neglected link 

between life and landscape. Geomorphology 157-158:6-16. 
Ward, J. V. 1994. Ecology of alpine streams. Freshw. Biol. 

32:277-294. 

Wallace, R. L., T. W. Snell, C. Ricci, and T. Nogrady. 2006. 

H. Segers, ed. Rotifera: volume 1: biology, ecology and 

systematics. 2nd ed. Backhuys Publishers, Leiden. 299 pp. 
Walther, G.-R., S. Beifiner, and C. Burga. 2005. Trends in the 

upward shift of alpine plants. J. Veg. Sci. 16:541-548. 
Webster-Brown, J., M. Gall, J. Gibson, S. Wood, and I. Hawes. 

2010. The biogeochemistry of meltwater habitats in the 

Darwin Glacier region (80 degrees S), Victoria Land, 

Antarctica. Antarct. Sci. 22:646-661. 
WGMS - World glacier monitoring services. 2008. WGMS 

Service, ed. Global glacier changes: facts and figures. 

Available at http://www.grid.unep.ch/glaciers (accessed 13 

September 2012). 
Woodward, G., D. M. Perkins, and L. E. Brown. 2010. Climate 

change and freshwater ecosystems: impacts across multiple 

levels of organization. Philos. Trans. R. Soc. Lond. B Biol. 

Sci. 365:2093-2106. 
Yeates, G. W., T. Bongers, R. G. M. De Goede, D. W. 

Freckman, and S. S. Georgieva. 1993. Feeding habits in soil 

nematode families and genera - an outline for soil 

ecologists. J. Nematol. 25:77-92. 

Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Table SI. Methods of preservation and laboratory treat- 
ment. 

Table S2. Taxon list. 

Table S3. PERMANOVA (pairwise) test results of 
seasonal differences in the MMB groups. 



2844 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



